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Abstract 

A new simple method for the first order phase transition kinet- 
ics is suggested. The metastable phase consumption can be imag- 
ined in frames of the modisperse approximation for the distribution 
of the droplets sizes. In all situations of the metastable phase decay 
this approximation leads to negligible errors in the total number of 
droplets appeared in the system. An evident advantage of the pre- 
sented method is the possibility to investigate the situation of the 
metastable phase decay on several sorts of heterogeneous centers. 

Attempts to give a theoretical description for the first order phase tran- 
sition appeared practically simultaneously with creation of the classical nu- 
cleation theory [|l| . An idea to give the global description of the first order 
phase transition was attractive and induced many publications (for example, 
see Wakeshima 0, Segal' Q). But all of them were intended only to estimate 
the main characteristics of the phase transition. The time of the cut-off of 
the nucleation clearly observed in experiments was adopted in these publica- 
tions without any proper justification. The first theoretical description with 
explicit determination of time behavior of supersaturation was given in 
where a homogeneous condensation was investigated. The method used in 
[|] was an iteration solution of integral equation of a substance balance. 

To give a global picture of phase transition one has to take into account 
a presence of active centers of condensation. The iteration method can be 
spread on heterogeneous condensation on similar centers (see [^), but for con- 
densation on several types of centers one can not calculate iterations with 



a proper accuracy (see also ||^). The system of condensation equations be- 
comes so complex that it can not be directly solved without simplifications. 
As the result it would be rather attractive to suggest some simple approxi- 
mations which can be used in the complex situation of the condensation on 
several types of centers. Certainly, this approximation has to be well based. 

Here we shall suggest a monodisperse approximation of the droplets size 
distribution to calculate the number of molecules in the liquid phase. From 
the first point of view this approximation is strange - we have already at- 
tracted attention to the necessity to know the behavior of a supersaturation 
which determines the form of the droplets size spectrum. But here we are 
going to show that with an appropriate choice of the cut-off (which can be 
also justified) one can give an adequate description of a nucleation period. 

The monodisperse approximation presented here differs from the total 
monodisperse approximation used in |J7[ for description of the intensive con- 
sumption of the metastable phase and can not be applied during a nucleation 
period. Here we use a special modification of the mentioned approximation 
which allows to describe the period of nucleation. 

This publication is intended to give the main idea of the monodisperse 
approximation which would be clear for experimenters. That's why we start 
from situations of homogeneous condensation and heterogeneous condensa- 
tion on similar centers which can be solved even without this approximation. 
Some technical details are excluded to give the most clear version (for exam- 
ple, a complete version of limit situations and monodisperse approximation 
in the intermediate situation is described in 0, the overlapping of the inter- 
mediate and limit situations is analysed in [Q, the transition of the special 
monodisperse approximation to the total one is discussed in 0). 

We use the physical model of nucleation kinetics described in |^. It is 
rather standard but to avoid misunderstanding we shall consider as the 
base of references. 

1 Homogeneous nucleation 

The condensation kinetics equation for the number G of the molecules in the 
liquid phase can be written in a well known form Q], 

G{z) = f r dx{z - xf exp(-rG'(a;)) (1) 
Jo 
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where parameter / is the amphtude value of the droplets sizes distribution 
r is some fixed positive parameter. 

One can analyse behavior of subintegral function g defined by 

G{z)^ rg{z,x)dx (2) 
Jo 

as a function of a size p = z — x, z. This function has the sense of the 
distribution of the number of molecules in droplets over their sizes p. 

In some "moment" 2; ( or t{z)) it can be presented in the following form 

• When p > 2; it is equal to zero (there are no droplets with such a big 
size) 

• When p < it is also equal to zero (there aren't droplets with a negative 
size) 

• At the intermediate p it grows rather quickly with a growth of p. It is 
easy to note that it grows faster than p^ grows. 

Really, if one takes into account that supersaturation decreases in time 
then we get g ~ p^. But supersaturation falls in time and there aren't so 
many droplets of the small size as of the big size. 

As the result one can see that the function g as the function of p has the 
sharp peak near p ^ z. This property takes place under the arbitrary z (or 
t{z)). 

The sharp peak of g allows to use for g the monodisperse approximation 
- a representation in the 6-like form with a coordinate corresponding to a 
position of the peak of function g, i.e. 

g ^5{p^ z) 

As the result one can state that the monodisperse approximation is based 
now. But it is necessary to determine the number of droplets in this approx- 
imation. 

It would be wrong to include the total number of already appeared 
droplets in this peak. Really, in the spectrum of sizes there are many droplets 
with small sizes. One can not describe these droplets as containing the same 
substance as the droplets of a big size. It would be more correct to exclude 
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them from the substance balance. So, it is necessary to cut off the small 
droplets. It can be done according to two recipes. 

The first recipe is the differential one. One can note that during all times 
which don't exceed the time of nucleation essentially the function g near 
maximum is close to 

Qappr f P 

This approximation corresponds to the constant value of supersaturation. 

One can cut off this approximation at a half of amplitude value (i.e. at a 
level fz^/2). Then one can get for the width Adif/z the following expression 

Adiffz = (1 - 2-V3)^ 

This cut off means that all droplets p < z — A^iffZ are excluded from 
consideration and all droplets with p > z — A^iffZ are taken into account in 
a (5-like peak. 

The second recipe is the integral one. One can integrate Qappr and require 
that ^ 

/ dappri^^i x)dx — N Z 
J 

An integration gives 

rz rz 

/ 9appr{z, X)dx = f {Z- xfdx = f 

Jo Jo 



z^ 



The width of spectrum is defined from condition that the number of droplets 
has to be equal to the amplitude multiplied by the width of spectrum Aintz: 

N = fAintZ 

This gives the following expression 

AintZ = z/4: 

One can see that Adif/z and AmtZ practically coincide. This shows the 
high selfconsistency of this approximation. The second recipe will be more 
convenient for concrete calculations. 

In fig.l one can see the application of the monodisperse approximation in 
the homogeneous case. 
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As the result one can say that all parameters of approximation are defined. 
Now it will be used to solve (|ID. 
Instead of (|I|) one can get 

G{z) = N{z/A)z^ 

where 

N{z/A) = fz/A 
is the number of droplets formed until t{z/A). This leads to 

G{z) = fzVA 

which coincides with the resulting iteration in the iteration method [^]. 
It is known (see 0,0) that this expression is very accurate which shows the 
effectiveness of the monodisperse approximation. Here the cut off of the tail 
of the sizes spectrum compensates the unsymmetry of the initial spectrum. 

The main result of the nucleation process is the total number of the 
droplets which can be found as 

Ntot = f dx exp{—TG{x))dx 
Jo 

Ntot = f / dxexpi-fTz^A) = f/^T-'/^D 
Jo 

where _ 

/•oo 

D = / exp(-xV4)rfx = 1.28 

JO 

The error of this expression is less than two percents (it is the same as in 
the iteration method). 

2 Heterogeneous condensation on similar cen- 
ters 

The condensation equations system can be written in the following form |^ 

G(z) = f r dx(z - X? exp(-TG(x))e(x) 
Jo 
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9{z) = exp(— 6 / exp{—TG{x))dx) 







with positive parameters /, b, T. An appearance of a new function 6 which is 
a relative number of free heterogeneous centers requires the second equation. 

The first equation of the system is rather analogous to the homogeneous 
case. The subintegral function here is also sharp. A function 6' is a decreasing 
function of time according to the second equation of the system. Then the 
function g which is again determined by (0) is more sharp than in the homo- 
geneous case. As far as the supersaturation has to fall one can see that g is 
more sharp than gappr- It allows to use here the monodisperse approximation 
for all z or t{z). 

As the result the monodisperse approximation is based for heterogeneous 
condensation. One needs here only the sharp peak of g{p) which can be easily 
seen. 

The successive application of the monodisperse approximation in the ho- 
mogeneous case shows that all droplets necessary for a metastable phase 
consumption at t{z) were formed until t{z/4:). In the heterogeneous case the 
exhaustion of heterogeneous centers increases in time. So, all essential at 
t{z) droplets were formed before t{z/4). 

At the same time the presence of a long tail in the situation of a weak 
exhaustion of heterogeneous centers requires to cut off the spectrum for the 
monodisperse approximation. As far as the long tail is essential in the situa- 
tion of a weak exhaustion one has to cut off the spectrum by the same recipe 
as in the situation of the homogeneous condensation: one has to exclude all 
droplets formed after z/4 which have the sizes p < z — AintZ = z — z/4. 

One can see in fig. 2 the monodisperse approximation in the situation of 
the heterogeneous condensation on similar centers. The form of the spectrum 
in this situation is illustrated in fig. 3. So, the way to construct approxima- 
tion is known. Now one can turn to concrete calculations. 

The number of the droplets formed until t{z/4) has to be calculated as 




An approximation for G has the form 



Giz) = L(l-9iz/A))z' 
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The total number of droplets can be determined as 



or 



iVtoi = ^(1-^(00)) 

f 

Jo 







or 



Ntot = ((1 - exp(-6^"exp(-r^(l - e{z/A))z')dz)) 



or 



poo 



rz/'l 



N^^^^ Ui-exp{-b exp(-rMl-exp(-6 / exp{-rG{x))dx))z^)dz)) 



The last expression has a rather complicate form. It contains several itera- 
tions in a hidden form which ensures the high accuracy. 

The last expression can be simplified. One of the possible recipes is the 
following. One can note that an expression for G is necessary at TG ~ 1. 
Then z attains some values A^^. But until ^(z/A the value TG is small and 
exp(rG(2;)) is close to unity. This leads to simplification of last expression 
which can be written in the following form 



Then one can fulfil calculation according to the last formula. The relative 
error is less than two percents. Here it is a little bit greater than in the 
homogeneous case because the form of initial spectrum is changed and there 
is no full compensation of the unsymmetry of spectrum and an exclusion of 
the tail. The relative error in the situation of heterogeneous condensation on 
similar centers is drawn in fig. 4. 

Now one can turn to explicit calculation of the integral in the last expres- 
sion. After the appropriate renormalization the subintegral function is more 
sharp than exp(— x^) and more smooth than exp{—x^). Both these functions 
have a sharp back front of spectrum. It allows to introduce the characteristic 
scale Az by equation 




|r(l - exp{-b{Az)/4:)){Azf ^ 1 
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Then 



Jo 



'OO 



exp(-r|(l 



exp{-bz/4:))z^)dz = Az 



A + B 
2 



where 



'OO 



A 



f 

Jo 



exp{-x^)dx = 0.89 



/' 

Jo 



oo 



B 



exp{—x'^)dx = 0.90 



Now the calculation is reduced to some algebraic manipulations. The 
error of the last approximation is less than one percent. 
As the result 



One can note that it is possible to formulate the recipe already in terms 
of Az. The long way is adopted here to give the most clear picture for the 
monodisperse approximation. 

3 Nucleation on several types of heteroge- 
neous centers 

The main advantage of monodisperse approximation is the possibility to use 
it for the condensation on the several types of centers. The iteration proce- 
dure can not be applied in this case successfully. The result of calculations 
according to shows this fact explicitly. The reason is the existence of the 
cross influence of the different types of centers through vapor consumption. 

In the condensation on similar heterogeneous centers in the situation of 
exhaustion the influence of this phenomena on the vapor consumption isn't 
important because in the situation of consumption the converging force of the 
heterogeneous centers exhaustion is extremely high. But in the situation with 
two types of heterogeneous centers the exhaustion of the first type centers 
can have a certain influence on a vapor consumption but the exhaustion of 
the second type centers is weak and there is no converging force due to the 
weak exhaustion of the second type centers. 

This effect in very thin and it can not be taken into account in the second 
iteration. But one can not calculate the third iteration analytically and this 
stops an application of iterations. Really, this phenomena isn't evident from 
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the first point of view but it exits and leads to tlie error of the second iteration 
in many times. 

The apphcation of the monodisperse approximation is based on the sharp- 
ness of function g. This property takes place already in this situation. So, 
there are no objections to apply the monodisperse approximation here. 

Here we shall reproduce the same formulas but with the lower indexes 
which determine the sort of heterogeneous centers. 

The system of condensation equations can be written in the following 



where the lower indexes denote the sorts of centers. This system can be seen 
by the direct generalization of the one type case. 

The subintegral function in the substance balance equations is also sharp. 
As far as all 9i are the decreasing functions of arguments then the function g 
defined by (0) (with proper indexes) is sharper than without the exhaustion 
of heterogeneous centers. So, due to the supersaturationdecreasing g is more 
sharp than gappr- It allows here to use the monodisperse approximation for 
all z or t{z). 

As the result one can see that the monodisperse approximation in this 
case is justified on the base of the sharpness of g{p). 

The same properties as in the previous case can be also seen here. One 
has to cut off the spectrum at z/A. Here all justifications are absolutely same 
as in the previous section. The characteristic situation for the nucleation on 
two types of heterogeneous centers is drawn in fig. 5. As the result the way to 
construct the monodisperse approximation is known. Now one can present 
calculations. 

The number of the droplets formed until t{z/4) on the centers of sort i 
has to be calculated as 



form 1^ 




Ar,(./4) = |(1- 0.(^/4)) 



An approximation for Gi can be now presented as 



G,{z) = ^-^{l-e,{z/A))z' 
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The total number of droplets is defined as 

Ni tot ^h^- Hoc)) 

Or 



or 



or 



or 



Ni tot = 7^(1 - exp(-6i / exp(-r ^^^(a;))^^)) 
Oi Jo j 

Ni tot = hi - eM-h /°°exp(-r5: ^(1 - 9j{z/A))z')dz)) 

Oi Jo , Oo 



j "J 



Ni tot = 7:^(1 - exp(-6i / exp(-r 
Oi Jo 

^^(1 -exp(-6, r^\xp{-TY,Gk{x))dx))z')dz)) 

Now one can simplify the last expression by the same way as in the one type 
case. 

Expressions for Gi arc essential at TJ2jGj ~ 1. Then z is near A^z. 
Until A(^z/4: the value T J2j Gj is small and exp(r J2j Gj{z)) is near unity. It 
leads to 

f. fOO f. 

Ni tot = r(l - exp{-bi / exp(-r^ (^(1 - exp(-6,-^/4))^3)^^)) 

Oi Jo J Oj 

Now one can fulfil the calculations according the explicit formula. The 
relative error of the last expression is less than five percents (here it increases 
slightly due to the complex form of the spectrums on different sorts. The 
relative error in the number of droplets is drawn in fig. 6. The calculation of 
the last integral is absolutely analogous to the previous section. The subin- 
tegral function after renormalization lies between exp(— a;^) and exp(— a;^). 
It allows to get the characteristic size Az from 

j 

Then 

exp(-r^ 1^(1 - eM-bjz/A))z')dz = ^z^^ 
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The relative error of the last expression is less than one percent. 
As the result 



The formula is similar to the final expression in the previous section. But 
parameters in the last formula have to be determined in another manner. 

The physical sense of the last expression is the separate exhaustion of 
heterogeneous centers. One sort of centers can influence on the other sort 
only through a vapor consumption. This fact can be seen also in the initial 
precise system of the condensation equations. 
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Fig.l 

Monodisperse approximation in homogeneous condensation. Here one can 
see four pictures for different periods of time (or for different values of z . One 
can introduce according to VG{/S.z) — 1 and it will he the characteristic 
scale of the supers aturation fall. In part "A" z = Az/2, in part "B" z = Az, 
in part "C" z = 3Az/2, in part "D" z = 2Az. One can see that the spectrums 
in part "A" and part "B" are practically the sam,e. It corresponds to the 
property of the similarity of spectrums until the end of the nucleation period. 
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Fig. 2 

Monodisperse approximation in condensation on the similar centers. The 
value Az is the same as in the previous figure ( i. e. it is determined without 
the exhaustion of centers). Now rG{Az) < 1 and Az will be the character- 
istic scale of the supers aturation fall in the situation without exhaustion. In 
part "A" z = Az/2, in part "B" z = Az, in part "C" z = 3Az/2, in part 
"D" z = 2Az. One can see that the spectrvms in part "A" and part "B" 
aren't similar. Now all spectrums are more sharp than in the homogeneous 
case. 
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Fig.3 

One can see two curves which are going from p — z to the small sizes 

p. The lower curve corresponds to the real spectrum calculated with account 
of heterogeneous centers exhaustion. The upper curve corresponds to the 
condensation without exhaustion of heterogeneous centers which is the worst 
situation where there is converging force due to the centers exhaustion. 

Concrete situation drawn here corresponds to b — 2 after renormalization 
(the values of parameters f and T can be canceled by appropriate renormal- 
ization). The value of z here equals to 3Az/2. 

The solid lines correspond to the precise numerical solution. The dashed 
lines correspond to application of monodisperse approximation. One can not 
separate the numerical solutions from the approximate ones except the slight 
deviation in the region of small p. As far as all (precise and approximate) 
solutions will go to zero there will be no deviations for z ^ Az (i.e. we stop 
at the worst moment). 
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o.oi: 



1 



Fzg.4 

The relative error of approximate solution for the nucleation on the sim- 
ilar heterogeneous centers. Here the values f and T can be canceled after 
renormalization and there remains only one parameter b which is the argu- 
ment of the function drawn here. It is clear that this function is small. All 
asymptotes can be checked analytically (see j^). 
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Fig.5 

Characteristic behavior of size spectrums for nucleation on two types of 
heterogeneous centers. One can cancel fi, T by renormalization. One can 
put fi < 1 due to the choice of centers. Here hi — 2, fi — 1/2, 62 = 1/2. 
The value z is taken as 2A2; (see fig.l). The are three curves here. The 
lower one corresponds to the spectrum of droplets formed on the first type 
centers, the intermediate one corresponds to the droplets size spectrum for 
the second type centers, the upper one corresponds to the spectrum calculated 
without exhaustion of heterogeneous centers ( the reasons are the same as in 
fig. 3). The solid lines are the numerical solutions, the dashed lines are the 
approximate solutions. One can not see the difference for the lower curve. 
For the upper and intermediate curve one can see only very slight difference. 
All spectrums are renormalized to have one and the same amplitude (which 
is marked by f). 
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Fig. 6 

Relative error for the nucleation on two types of centers. For two types of 
centers there exists five parameters (two of them can he canceled by renormal- 
ization). We have already adopt that T is one and the same for different types 
of centers (for the reasons see j^, for numerical results see j^, for analytical 
estimates see /|T^/, for recipes of calculation in this situation see /[H|// We 
cancel here fi, T. We consider /2 < 1 and the first component is the leading 
component in the metastatic phase consumption. So, the worst situation for 
the error in the droplets number formed on the second type centers (which is 
drawn here ) will be when 62 = and the is no converging force of the centers 
exhaustion. So, there remain two parameters f2 (it is marked by f) and bi 
(which is marked by b). One can see that the error is small. The calculations 
for the variations of all parameters give the same value of error but it is hard 
to reproduce these numerical results. 
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